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Abstract. The lattice Gel'fand-Dikii hierarchy was introduced in ) 131 as the 
family of partial difference equations generalizing to higher rank the lattice 
Korteweg-de Vries systems, and includes in particular the lattice Boussinesq 
system. We present a Lagrangian for the generic member of the lattice Gel'fand- 
Dikii hierarchy, and show that it can be considered as a Lagrangian 2-form when 
embedded in a higher dimensional lattice, obeying a closure relation. Thus the 
multiform structure proposed in [8] is extended to a multi-component system. 



1. Introduction 

Multidimensional consistency [HJ E] has come to be regarded as one of the hallmarks 
of integrability for discrete systems. In brief, it is the property that several copies of 
an equation may be imposed simultaneously on a higher dimensional lattice, and no 
inconsistency or multivaluedness occurs in the evaluation of the dependent variables 
at each lattice site. Lagrangian multiforms were proposed for multidimensionally 
consistent discrete systems in [S] to remedy the fact that while several equations 
may coexist on the multidimensional lattice, only one equation can be derived from 
a scalar Lagrangian. The multiform structure allows copies of the relevant equation 
in all possible lattice directions to be derived through the Euler-Lagrange equations, 
making multidimensional consistency manifest on the level of the Lagrangian. 

The lattice systems considered in the first instance were quadrilateral equations 
of the form 

Q{u,u ll ,'u v ,Uf l v)a lll ,a v ) = Q, (1.1) 



where u — u{n^n v ) depends on two discrete variables n^,n v , shifts of u in the n 
direction are denoted by it M (so that for example u M = u(n^ + l,n v )), and the 
are lattice parameters associated with the n^-direction. The examples, defined on 
elementary plaquettes in a 2-dimensional lattice, were chosen from the classification 
given in [T] . 

The important observation was that all of the systems considered in [8] admit 
Lagrangians = L^yiu^^u^a^au) in terms of 3 points, which satisfy the 
following closure relation 

A x £»u + + A„£a m = 0, (1.2) 

where the difference operator A a acts on functions / of u = u(nA,n M ,n y ) by the 
formula A\f(u) = f(u\) — f(u), and on a function g of u and its shifts by the formula 
A\g(u, u^, u v ) — g{u\, itAju, u\ v ) — g(u, u^, u w ). This allowed the interpretation of the 
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Lagrangian as a closed 2-form on the multidimensional lattice, and on the basis of this 
a new variational principle was proposed. 

This idea was extended in [§] to a 3-dimensional system, the bilinear discrete 
Kadomtsev-Petviashvili (KP) equation. This is an equation in 6 points, with a 
Lagrangian C\^ v also in terms of 6 points. In this case the relevant closure relation 
obeyed by the Lagrangian is 

A p £a^j/ — &-\£-nv P + Ap£„ P A — A U C P \ P = 0, (1.3) 

where as before the difference operator A\ acts on functions / of r = t(ti\, n p ,n v , n p ) 
by the formula A\f(r) — f{r\) — /(r), and on a function g of r and its shifts by the 
formula A\g(r, r M , r v , t p ) = g(r x t\ p ) - g(r, , Ty , Tp). Here, the Lagrangian 

can be considered as a closed 3-form. 

So far, all systems with Lagrangians that have been shown to obey a closure 
relation have involved variables around a simple plaquette or cube. This is also the case 
for the lowest order member in what is effectively a lattice analogue of the Gel'fand- 
Dikii (GD) hierarchy [T3], the discrete Korteweg-de Vries (KdV) equation, which was 
one of the systems considered in |S]. The next system in the hierarchy is the discrete 
Boussinesq equation [13] . this involves not only nearest neighbouring points, but also 
next-nearest neighbours. Along with higher order systems in the GD hierarchy, it can 
be written as a coupled system of partial difference equations defined on an elementary 
plaquette. In this paper we present a Lagrangian for the generic member of the lattice 
GD hierarchy, and show that it obeys the closure relation (| 1 . 2|) , allowing us to apply 
the variational principle proposed in [8]. 

2. The lattice Gel'fand-Dikii hierarchy 

The lattice GD hierarchy first appeared in [13] , where the direct linearization method 
was used to find a discrete analogue of the continuous GD hierarchy, which is a 
hierarchy of systems associated with higher order spectral problems [101 |4j . As 
already mentioned, the first members in the hierarchy are the lattice KdV and lattice 
Boussinesq equations, higher order members are coupled systems of partial difference 
equations in terms of variables u, Vj, Wj, where < j < N — 2, given by the following. 

Vj+l — Vj+i = (p — q + u — u)vj — pvj + qvj (2.4a) 

Wj+i ~ Wj+x = — (p — q + u — u)uij — qu)j + pulj (2.4b) 
for < j < N - 3, and 

(jp-q + u- u)(v N -2 - w N - 2 ) = (p + q + u)[(p - q + u - u)v N ~ 3 -pv N ^ 3 + qv N - 3 ] 

N-3 

- Wj ((-p) N -^u-(-q) N -^u)] 
N-2 N-l-j 

"EE WiU-P)"- 1 -^-!-^)"- 1 -^-!], 
3=2 t=0 

(2.4c) 
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identifying v — Wq 



Here we have used the notation we find most instructive: 



the dependent variables are u,Vj,Wj for < j < N — 2, and we consider them to 
depend on two independent variables n, m. The symbol ~ is used to denote shifts in 
the n-direction, and ^ denotes shifts in the m-direction, so that if u = u(n, m), then 
u = u(n + 1, m) and u = u(n, m + 1). The lattice parameters p, q are associated with 
the n, m-directions respectively. This is illustrated in the diagram below. 
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The variables Vj, uij are evaluated at the same lattice point as u, while Vj, Wj are 
evaluated at the same lattice point as u, and so on. 

As noted in [13] the lattice GD hierarchy arises from a Zakharov-Shabat type of 
linear problem 



in which 



(p + uk)(p k = L k ■ ip k , (q + Lok)(p k = M k ■ (p k 

\ 

-vi p 1 

L k = 



VN-2 
N 



o 

* WN-2 



P 1 
Wi p + u J 



(2.5) 



(2.6) 



and where the matrix M k is a similar matrix obtained after the replacements p i— ► q 
and ~ i— »^ . The term * in the left-lower corner of the matrix L k is such that the 

x (— k) , i.e. we have the expression 



determinant &ei{L k 

N-2 



J2 (-p) n -^ ® 

3=0 



N-l N-l-j 



WiVj-i 



(2.7) 



Because the system is multidimensionally consistent, a concept introduced 
independently in [14] and [2], we are free to impose copies of the equations in other 
lattice directions, with appropriate lattice parameters. In particular, suppose we have 
another lattice direction associated with parameter r, where shifts in this direction 
are denoted by ~ . Then copies of the equation (|2.4aj) can be imposed on each pair of 
lattice directions, giving 



Vj+l - Vj+i = (p - q + u — u)Vj — pdj + qvj 
Vj + i — Vj+i = (q — r + u — u)Uj — qUj + rvj 
Vj+i — Vj+i = (r — p + u — u)vj — rvj + pUj 



(2.8a) 
(2.8b) 
(2.8c) 
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for < j < N - 3. 

In particular, summing equations (|2.8ap . (|2.8bp and (|2.8c|l for ,7 = gives 

= (p — q + u — u)(u + r) + (q — r + u — u)(u + p) + (r — p + u — u)(u + q), (2.9) 

which is in fact the lattice Kadomtsev-Petviashvili (or lattice KP) equation [11]. That 
(|2.9p holds is natural, since members of the lattice GD hierarchy can be viewed as a 
special type of periodic reduction of the lattice KP equation. 



3. Lagrangian for the lattice GD hierarchy 

Lagrangians for the two lowest order members of the lattice GD hierarchy have already 
appeared in the literature, an action for the KdV lattice was first given in [3], and for 
the Boussinesq lattice in [T3] (these systems arise from the lattice GD hierarchy by 
taking N — 2 or N — 3 respectively). In fact, as we show here, it is possible to write 
a Lagrangian for the generic member of the hierarchy. 
Note first that if we define 

7j (p, q) EE ( ~ P)J+1 ~ }~ q)3+1 = (-1)^ +P , -\ + • • ■+pq j - 1 + q J ), (3.10) 

-p + q 

then equation l|2.4c[) can be written in a more convenient form as 

/ \ / \ N-3N-3-i 

, s . ~ (p- q)iN-i{p-,q) ~ 

lN-l(P,q) + VN-2 - WN-2 = ~ h > > 7N-3-i-jVjWi 

p — q + u — u 

y H t=0 j=0 

N-3 

~ 7N-2-j (Vj ~ Wjj , (3.11) 

3=0 

which enables us more easily to see that the following proposition holds. 

Proposition 1: The system consisting of equations (I2.4ap . (|2.4b|) and (|3.1ip solves 
the discrete Euler- Lagrange equations for the following Lagrangian 

£pq = (p-<l)lN-i(p,q)ln(p-q + u-u) -j N -i(p,q)(u-u) 
- lN-2-j(p,q)(u " u)Vj 

3=0 
AT-3JV-2-J 

- 22 X/ lN-2-i-j (p, q)wi [vj -Vj~(p-q + u- u)vj-i + pvj-x - qvj-x] 

i=0 j=l 

(3.12) 

under independent variation of u,v and w. 

Proof: Here we take the usual point of view and consider the action to be the 
sum of the Lagrangians over all n, m, i.e. 
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The discrete Euler-Lagrange equations arise as a consequence of the requirement that 
SS = 0. We have 

o = ss 

= V i — iN-i(p,q) {ou - <5u) 

^ I p-q+u-u 

AT-2 
AT-3 N-2-i 

- iN-2-i-j (p, q){Swi [vj -Vj - (p-q + u- u)vj-i + pvj-i - qvj-i] 

i=0 j=l 

+ u>i[<5u;' — (JtJj — (<5u — Su)Vj—\ — (p — q + u — — q8vj-i\} 

_ y-' / / (p - q)iN-i(p,q) (p- q)iN-i(p,q) 

„ez^ p-q + u-u p-q + u-u 

N-3N-3-i ^ _ 

+ X! X! jN-3-i-j(p,q)(v j W l -VjWi) 
i=0 j=0 

N-3 ^ ^ 

- W-2-j(p, q){vj - vj - wj + wj) 

- V N _ 2 + WAT-2 + WN-2 - WN~2^J Su 
N-3 N-3-i 

~ lN-3-i-j (P, q) fij+l - V j+ i - (p-q + u- u)Vj + pVj - qv J ]Sw t 

i=a j=o 

N-3N-3-i 

+ Y ^ lN-3-i- 3 (p, q) [w%+l ~ Wi+l + (p- q + u- u)w l + qw t - pWi]§Vj > . 
i=a j=o ' 

(3.14) 

The coefficient of Swi must be zero, which gives us equation (|2.4ap , the coefficient 
of Svj must be zero, which gives us equation (|2.4b[) . and the remaining term which 
multiplies Su is two shifted copies of equation p. lip . Hence the lattice Gel'fand-Dikii 
system of equations solve the Euler-Lagrange equations for the Lagrangian (|3.12p . ■ 



At this point we would like to make several remarks. 

(i) For clarity we have chosen to label the Lagrangian with the lattice parameters p, q 
to indicate it is defined on a plaquette in a 2-dimensional surface corresponding 
to the respective lattice directions. 

(ii) The Lagrangian (|3. 1 2[) is antisymmetric with respect to the interchange of the 
lattice directions associated with the parameters p, q, a property which will be 
important when we come to define the multiform structure. 

(iii) Total derivative terms have been included which at first sight may appear 
superfluous. They do, however, prove necessary in the verification of the closure 
relation below. For the same reason, we do not have the freedom to multiply 
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the Lagrangian by constants involving the lattice parameters p, q, we may only 
multiply by true constants. 

The main result of the paper is the following. 

Proposition 2: The Lagrangian defined by (|3.12p satisfies the following closure 
relation on solutions to the lattice GD hierarchy equations when embedded in a 3- 
dimensional lattice. 

ApCq r + AqC r p + A r Cpq — 0, (3.15) 

where the difference operator A r acts on functions f of u — u(n p ,n q ,n r ) by the 
formula A r f{u) — f(u) — f{u), and on a function g of u and its shifts by the formula 

A r g(u, u, u, u) — g(u, u,u,u) — g(u, u, u, u). 

Proof: Firstly, on equation (|2.4a|) is is clear that the last term in the Lagrangian 
will disappear. This leaves us with 

r = T + c + c — c - c — c 

— *^pq * *^qr <^ *^rp *~*pq ^qr *~*rp 

i \ / m (p — q+u-u\ (q-r + u-u 

= (P-Q nN-i{P,q)ln[ — — - + {q-r)jN-i(q,r)lii[ — — - 

\p — q + u — u J \q — r + u — u 

(v — p ~\~ u — TL \ ~ 
— - — = -7N-i(p,q)(u-u-u + u) 
r — p + u — u J 

— lN-x(<lt r )(u — u — u + u) — 7jv_i(r,p)(u — u — u + u) 

N-2 , ^ ^ 

\ 7JV-2-i(p, q)[{u - u)Vj -(u- u)vj] + 7jv-2-j(<7, r)[(u - u)v 3 -(u- ujVj] 



3=0 



+ 7W-2-J (r, p) {(u - u)v 3 - (u - u)vj}>, (3.16) 



which, on rearranging the terms, is 



r = - {-p) n In 



p — q + u — u\ I r — p + u — u 



p- q+u-uj Vr -p+u-u 



_ f q \N ln ( ( q-r + u-u \ ( p-q+u-u 
\\q - r + u -u ) \p - q + % 

, \N i f f r ~ P + u ~u\ / q — r + U — u\\ ^ ~ _ „ 

- \-r) In — - — - — s - -f N -i{p,q){u -u-u + u) 

\\r -p + u-u J \ q -r + u-u) J 

— 7N-i(q, r)(u — u-u + u) - jN-i(r,p)(u — u — u + u) — (u- u)vn-2 

+ (u — u)v N -2 - (u — u)vn-2 + (u- u)vn-2 - (u — u)vn-2 + (u — u)vn-2 



N-3 

\ lN-2- 3 (p, q)[(u - u)Vj - (u - u)vj] + 7JV-2-j(<?, r)[(u - u)vj - (u- 



+ lN-2- {r,p){{u - u)vj - (u - u)vj] >. (3.17) 



We have already shown in section 2 that the lattice KP equation (|2.9p holds provided 
that equation (|2.4a|) holds. Using this fact, it is clear that the logarithm terms 
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disappear. So we are left with 

F = — 7jv_i(p, q)(u — u — u + u) — jn-i(q, r)(u — u — u + u) — -fN-i{r,p)(u — u — u + u) 

+ u(^N-2 ~ VN-2) + u(vn-2 - V N - 2 ) + u(y N -2 - VN-2) 

W-3 , _ _ „ 

~ ^2 } TiV-2-j(P) q)[(u - u)Vj -(U- UjVj] + lN-2-j(q, r)[(u - u)Vj -(u- u)Uj] 

+ ^N-2-]{r,p)[(u-u)v ] - (u -¥)«,•] I . (3.18) 
Here it is helpful to introduce a new object 

( ] = ~ _ ^( 7 J'+ 1 (' ? ' r ) _ 7j+iO*,P)) ; ( 3 - 19 ) 

which is invariant under cyclic permutations of p, q, r. This is not immediately 
apparent, but is due to the fact that -fj obeys the relation 

{p - q)lj {P, q) + (q- r)jj (q, r) + (r- p)jj (r, p) = 0, (3.20) 

which allows us to write 



1 


p 




q 




1 




p 




q 




1 




r 




p 



7j+i(p,q)-7j+i(q,r)J, (3.21) 

and this is clearly (|3 . 1 0|1 after a cyclic permutation of p, q and r. Hence ej is invariant 
under such cyclic permutations. The following identity for ej also holds 



e i+1 + re,j = j j+1 (jp, q), (3.22) 

lj+i(q,r) +nj(q,r) 



since 

(- q y+ 2 - (- r y+ 2 (-qy +1 - (-ry +1 



-q + r -q + r 

= (-q) j+ \ (3.23) 
and similarly jj+i(r,p) + r~fj(p, q) = (— p) 3+ , so that 

e j+1 + re 3 = — !— ^7 i+2 (g, r) - 7 J+2 (r,p)^ + — ^— (jj+i{q, r) - 7j+x(r,i>)^ 

= -^—^((l3+2(q,r) + n J+1 {q 7 r)) - (j j+2 (r,p) + r-y j+1 {r,p)) 

_ (-q) 3+2 - i-p) 3+2 



p-q 

= 7j+i(P,?)- (3-24) 
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From the definition, = I, and so we can write 

F = — 7jv_i(p, <?)(u — u — u + u) — 7jv-i(<7, r)(u - u — u + u) - -fN-i{r,p)(u — u — u + u) 
+ e a u{v N -2 - v N - 2 ) + e u(vN~2 - «jv-a) + tou(v N - 2 - v N - 2 ) 

N-3 , _ ^ ^ 

- ^2 } TiV-2-i(P) ?)[0" - - (U - UjVj] + lN-2-j(q, r)[(u - u)Vj -(u- u)Vj] 

3=0 

+ jN-2-j (r, p) [(u - u)vj - (u - ufUj] I . (3.25) 
For any < k < N — 3, the expression 

Afe = £ k u(vN-2~k - VN-2-k) + tku(vN-2~k ~ VN-2-k) + e k u(VN-2-k ~ UjV-2-fc) 
N-3-k f _ _ „ 

- \~/N-2-j(p,q)[(u -u)vj - {u-ufuj] +j N - 2 -j(q, r)[(u-u)vj - (u-u)%] 

+ jN-2-j(r,p)[(u - u)Vj - (u-u)Vj]^ (3.26) 
can be written as 

A fe = e k u{{p -q + u- u)v N ^ 3 ^ k - pv N - 3 ^ k + qv N - 3 -k) 

+ e k u((q -r + u- w)tJjv-3-fc - qvN-3-k + rv N - 3 -k) 
+ e k u((r -p + u~ u)v N - 3 ^ k - rujv-3-fc + pv N ^ 3 - k ) 

- 7fe+i(p, q)[(u - u)v N -3-k - (u - 2)wjv-3-fc] 

- 7/s+i (g, r)[(u - u)v N - 3 - k - (« - w)wAr_ 3 _fc] 

- 7fc + i (r, p) [(u - u)v N - 3 - k - (u - H)v N - 3 - k ] 

JV-4-& /• ^ _ 

- 51 I 7JV-2-j(p,t?)[(u ~ w)«j - («- +7N-2-j(q,r)[(u-u)vj - (u-u)i7,-] 

+ jN~2-j(r,p)[(u -u)Uj - (u - tZ)Uj]j, (3.27) 

where we have made use of (|2.4al) to eliminate the terms involving shifts of u/v-2-fc- 
On rearranging, 

Afe = [e k u(p — q + u — u) + e k u(q — r + u — u) + e k u(r — p + u — u) 

- Jk+i(p, q)(u -u)- Jk+i(q, r)(u -u)- j k+1 (r,p)(u - u)]v N - 3 - k 
+ u[(-qe k +'j k+ i(r,p))v N - 3 -k - (-re k +7 fe+ i(p, q))v N - 3 - k ] 
+ u[(-re k +7fe + i(p, q))vif-3-k - (~Pe/s + lk+i(q, r) )v N - 3 - k ] 
+ u[(-pe k +"f k +i{q,r))v N - 3 - k - {-qe k + jk+i(r,p))vN-3-k] 

N-i-k , ^ ^ 

- ^2 \ lN-2-j{p-,q)[{u -u)v 3 ~ (u-ujvj] +7iv_2-j(g,r)[(u - u)vj - (u-u)%] 

+ jn-2-j (r, p)[{u- u)Vj - (u - u)vj] \ , (3.28) 
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and then we can use (|2.9p on the very top line to give 
A fe = [(re* - j k+1 (p,q))(u -u) + (pe k - j k+1 (q,r))(u-u) 
+ {pe k - 7fc+i(r,p))(u -u)]v N - 3 - k 

+ u[(-qe k + "fk+i{r,p))vN-3-k - {-reu + 7fc+i(p, q))vN-3-k] 

+ u[(-re k +j k+ i(p, g))ujv-3-fc - (~pe/s + 7*H-i(3 r > r ))^JV-3-fc] 

+ «[(— pe* +7/c+i(?,? , ))wAr-3-fc - (-qek +jk+i(r,p))vN-a-k] 

JV-4-fe j- ^ 

- X! I "fN-2-j(p,q)[(u -u)Vj - (u - +7jV_2-j(g,r)[(t4-u)Uj - (u - u)Vj] 
3=0 { - 

+ n fN-2-j(r,p)[(u - u)Vj - (3.29) 
Here the identity (|3.22p comes into play, to give 

Afe = £fe + itt(lJiV-3-fc - VN-3-k) + £k+lu(vN-3-k ~ UjV-3-fc) + efc+l"(wA'-3-fc - VjV-3-fc) 
JV-4-fc j- ^ _ 

- I 7JV-2-j(p,g)[(« - w)tJj - (u - uJUj] +"fN-2-]{q,r)[(u - u)Vj - (E- tt)%] 

+ jN-2-j{r,p)[(u - u)vj -(u- u)vj] 

= Afc+i, (3.30) 

which means that for any < j, k < N — 2 we have Aj = Afe . This allows us to greatly 
simplify T, since 

F = — 7at_ 1 (p, q)(u — u - u + u) - 7Ar_i(g, r)(u - u — u + u) - 7Ar_i(r,p)(u — u - u + u) 
+ A 

= — 7iV-l(p, g)(w - u - u + u) - jN-i{q, r)(u -u — u + u) - -fN-i(r,p)(u — u-u + u) 
+ 

= — 7jv_i(p, — u — 2 + u) — jN-i{q, r)(u -u — u + u) - "f N -i(r,p)(u — u-u + u) 
+ eN~2u(u -u)+ €n~2u(u — u)+ €n-2u(u — u). (3.31) 
Using once again the equations (|2.9p and then (|3.22p . this is 

r = - 7at_i(p, q)(u — u — u + u) — jN~i(q, r)(u -u — u + u) - 7Ar_i(r,p)(w — u-u + u) 
+ reN-2(u — u — u + u) + pejy_2(w — u — u + u) + qe^^iu — u — u + u) 
= (re N -2 -jN-i(p,q)){u-u-u + u) + (pe N - 2 - 7jv-i(g,r))(u - u-u + u) 
+ {qe-N-2 ~ jN-i(r,p))(u -u - u + u) 

= — 6N-l(u — U — U + U,) — 6N-l(u — U — U + u) — £N-l{u — U — U + u) 

= 0. (3.32) 
Thus the closure relation is verified. ■ 



Note that in the above computation, the only equations used were the equation 
involving the Vj (|2.4a|l . and the lattice KP equation (|2.9p . the latter was shown earlier 
to be a consequence of copies of (|2.4aj) . 
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4. Lattice Boussinesq and KdV equations 

The lattice Boussinesq equation deserves special mention as it has attracted much 
interest lately, for example with regard to the Pentagram map [15] . It is a particular 
case of the lattice GD hierarchy, taking N = 3, and as such can be written as a system 
of equations in the variables u, t>i , i£>i . However, it is possible to eliminate v\, w\ and 
express the equation in terms of the variable u only, as follows. 

3 3 3 3 ^ 

— ^| — ~ ^—Si — 2 + b + 2<?)(M + w) - (2p + q)(u + u) 

p — q + u — u p—q+u—u 

+ (p — q + u — u)u + (p — q + u — u)u + uu — uu — 0. (4.33) 

Starting from the action given in [13] . we need to make only minor modifications in 
order to arrive at a Lagrangian in terms of the variable u only which satisfies the 
closure relation (|3.15p on solutions to the lattice KP equation (|2.9|) . 

Cpq = (p 3 — q 3 ) ln(p — q + u — u) - (p 2 + pq + q 2 )(u - u) + (p + q)(u - u)u 

+ (p — q + u — u)uu — puu + quu. (4.34) 

That the above Lagrangian satisfies the closure relation (|3 . 1 5[) can easily be verified 
by direct computation. 

The lattice KdV equation, which is the member of the lattice GD hierarchy where 
N = 2, is a more degenerate case and needs to be treated separately. Here we do 
not have the equations (|2.4aj) and (|2.4bj) as N is too small, we have only the equation 
(|3~TTj) . which is 

(p + q + u — u) (p — q + u — u) = p 2 — q 2 . (4.35) 
A Lagrangian was first given in [3] , which is equivalent to the following Lagrangian 

C pq = —(p 2 — q 2 ) ln(p — q + u-u) + (u- u){p + q — u). (4.36) 

Again, the closure relation holds on solutions of the lattice KP equation (|2.9|) . but we 
need to use copies of equation (|4.35[) in 3 lattice directions to show that (|2.9[) does 
indeed hold. The lattice KdV equation is a case already treated in [8]. 

5. Lagrangian multiform structure 

Having established a Lagrangian for each member of the lattice GD hierarchy satisfying 
the closure relation (|3 . 1 5[) . we now interpret this result in terms of a Lagrangian 
multiform structure. In fact, the existence of the closure relation allows us to develop 
the variational principle proposed in [8] for this class of systems. This comprises the 
following. 

Noting that the Lagrangian (|3.12[) is defined on an elementary plaquette, we can 
define an action S for any given surface a consisting of a connected configuration 
of elementary plaquettes o~ pq in the multidimensional lattice (where the labelling by 
the lattice parameters p, q indicates a vq lives on the sublattice corresponding to the 
respective lattice directions) by summing the Lagrangian contributions from each of 
the plaquettes in the surface cr, i.e. 



S[u,Vi > ...,v N -2,Wi,...,w N -2]a] = \ C pq , 



(5.37) 
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taking into account the orientation of the plaquette (as noted earlier, C pq has the 
property of antisymmetry with respect to interchange of the two lattice directions, so 
this sum is well-defined). This action depends not only on the dependent variables, but 
also on the geometry of the independent variables. The variational principle proposed 
in [8] amounts to the following reasoning. Imposing independence of the action under 
local variations of the surface, keeping the boundary fixed, requires the closure relation 
to hold. Furthermore, the surface independence allows us to locally deform the surface 
in any way we choose away from the boundary. In particular, away from the boundary 
we may render it locally flat so that we have here a regular 2-dimensional lattice 
on which we can apply the variational principle leading to the usual discrete Euler- 
Lagrange equations. It then follows from the proof of Proposition 2 in the previous 
section that for this specific Lagrangian, the equations of motion ()2.4a|) . (j2.4bp and 
p. lip are compatible with the surface independence. We view this circular mechanism 
behind the variational principle as a manifestation of multidimensional consistency on 
the level of the Lagrangian. 

The requirement that the closure relation should hold specifies to some extent the 
Lagrangian. For all examples in [SJ, for the closure relation to hold it was necessary 
for the Lagrangian to be antisymmetric with respect to the interchange of the two 
lattice directions, although total derivatives could be added so long as they were 
also antisymmetric. Here, however, it seems we are even more constrained, as in 
general we are not free to add such total derivatives. That the closure relation places 
such restrictions on the Lagrangian could be of relevance in the inverse problem of 
Lagrangian mechanics, a field of study which dates back to the 1880s j7], cf. [18] for 
a review. 

It would also be interesting to see whether a similar closure relation holds for 
the continuous GD hierarchy. A Lagrangian for the generating partial differential 
equation (PDE) for the KdV hierarchy [12] was shown to obey a continuous analogue 
of the closure relation in [8] . The generating PDE for the Boussinesq hierarchy along 
with its Lagrangian appeared in |16[ 117] , it is expected that this would also obey a 
continuous analogue of the closure relation. 
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